Automatization of metabolite extraction for high-throughput metabolomics: case study on transgenic isoprene-emitting birch.

Metabolomics studies are becoming increasingly common for understanding how plant metabolism responds to changes in environmental conditions, genetic manipulations and treatments. Despite the recent advances in metabolomics workflow, the sample preparation process still limits the high-throughput analysis in large-scale studies. Here, we present a highly flexible robotic system that integrates liquid handling, sonication, centrifugation, solvent evaporation and sample transfer processed in 96-well plates to automatize the metabolite extraction from leaf samples. We transferred an established manual extraction protocol performed to a robotic system, and with this, we show the optimization steps required to improve reproducibility and obtain comparable results in terms of extraction efficiency and accuracy. We then tested the robotic system to analyze the metabolomes of wild-type and four transgenic silver birch (Betula pendula Roth) lines under unstressed conditions. Birch trees were engineered to overexpress the poplar (Populus × canescens) isoprene synthase and to emit various amounts of isoprene. By fitting the different isoprene emission capacities of the transgenic trees with their leaf metabolomes, we observed an isoprene-dependent upregulation of some flavonoids and other secondary metabolites as well as carbohydrates, amino acid and lipid metabolites. By contrast, the disaccharide sucrose was found to be strongly negatively correlated to isoprene emission. The presented study illustrates the power of integrating robotics to increase the sample throughput, reduce human errors and labor time, and to ensure a fully controlled, monitored and standardized sample preparation procedure. Due to its modular and flexible structure, the robotic system can be easily adapted to other extraction protocols for the analysis of various tissues or plant species to achieve high-throughput metabolomics in plant research.

Pipecolic acid synthesis is required for systemic acquired resistance and plant-to-plant-induced immunity in barley.

Defense responses in plants are based on complex biochemical processes. Systemic acquired resistance (SAR) helps to fight infections by (hemi-)biotrophic pathogens. One important signaling molecule in SAR is pipecolic acid (Pip), accumulation of which is dependent on the aminotransferase ALD1 in Arabidopsis. While exogenous Pip primes defense responses in the monocotyledonous cereal crop barley (Hordeum vulgare), it is currently unclear if endogenous Pip plays a role in disease resistance in monocots. Here, we generated barley ald1 mutants using CRISPR/Cas9, and assessed their capacity to mount SAR. Endogenous Pip levels were reduced after infection of the ald1 mutant, and this altered systemic defense against the fungus Blumeria graminis f. sp. hordei. Furthermore, Hvald1 plants did not emit nonanal, one of the key volatile compounds that are normally emitted by barley plants after the activation of SAR. This resulted in the inability of neighboring plants to perceive and/or respond to airborne cues and prepare for an upcoming infection, although HvALD1 was not required in the receiver plants to mediate the response. Our results highlight the crucial role of endogenous HvALD1 and Pip for SAR, and associate Pip, in particular together with nonanal, with plant-to-plant defense propagation in the monocot crop barley.

Isoprene-Emitting Tobacco Plants Are Less Affected by Moderate Water Deficit under Future Climate Change Scenario and Show Adjustments of Stress-Related Proteins in Actual Climate.

Isoprene-emitting plants are better protected against thermal and oxidative stresses, which is a desirable trait in a climate-changing (drier and warmer) world. Here we compared the ecophysiological performances of transgenic isoprene-emitting and wild-type non-emitting tobacco plants during water stress and after re-watering in actual environmental conditions (400 ppm of CO2 and 28 °C of average daily temperature) and in a future climate scenario (600 ppm of CO2 and 32 °C of average daily temperature). Furthermore, we intended to complement the present knowledge on the mechanisms involved in isoprene-induced resistance to water deficit stress by examining the proteome of transgenic isoprene-emitting and wild-type non-emitting tobacco plants during water stress and after re-watering in actual climate. Isoprene emitters maintained higher photosynthesis and electron transport rates under moderate stress in future climate conditions. However, physiological resistance to water stress in the isoprene-emitting plants was not as marked as expected in actual climate conditions, perhaps because the stress developed rapidly. In actual climate, isoprene emission capacity affected the tobacco proteomic profile, in particular by upregulating proteins associated with stress protection. Our results strengthen the hypothesis that isoprene biosynthesis is related to metabolic changes at the gene and protein levels involved in the activation of general stress defensive mechanisms of plants.

Predicting functions of putative fungal sesquiterpene synthase genes based on multiomics data analysis.

Sesquiterpenes (STs) are secondary metabolites, which mediate biotic interactions between different organisms. Predicting the species-specific ST repertoires can contribute to deciphering the language of communication between organisms of the same or different species. High biochemical plasticity and catalytic promiscuity of sesquiterpene synthases (STSs), however, challenge the homology-based prediction of the STS functions. Using integrated analyses of genomic, transcriptomic, volatilomic, and metabolomic data, we predict product profiles for 116 out of 146 putative STS genes identified in the genomes of 30 fungal species from different trophic groups. Our prediction method is based on the observation that STSs encoded by genes closely related phylogenetically are likely to share the initial enzymatic reactions of the ST biosynthesis pathways and, therefore, produce STs via the same reaction route. The classification by reaction routes allows to assign STs known to be emitted by a particular species to the putative STS genes from this species. Gene expression information helps to further specify these ST-to-STS assignments. Validation of the computational predictions of the STS functions using both in silico and experimental approaches shows that integrated multiomic analyses are able to correctly link cyclic STs of non-cadalane type to genes. In the process of the experimental validation, we characterized catalytic properties of several putative STS genes from the mycorrhizal fungus Laccaria bicolor. We show that the STSs encoded by the L.bicolor mycorrhiza-induced genes emit either nerolidol or α-cuprenene and α-cuparene, and discuss the possible roles of these STs in the mycorrhiza formation.

Interaction with the entomopathogenic fungus Beauveria bassiana influences tomato phenome and promotes resistance to Botrytis cinerea infection.

Plants are central to complex networks of multitrophic interactions. Increasing evidence suggests that beneficial microorganisms (BMs) may be used as plant biostimulants and pest biocontrol agents. We investigated whether tomato (Solanum lycopersicum) plants are thoroughly colonized by the endophytic and entomopathogenic fungus Beauveria bassiana, and how such colonization affects physiological parameters and the phenotype of plants grown under unstressed conditions or exposed to the pathogenic fungus Botrytis cinerea. As a positive control, a strain of the well-known biocontrol agent and growth inducer Trichoderma afroharzianum was used. As multitrophic interactions are often driven by (or have consequences on) volatile organic compounds (VOCs) released by plants constitutively or after induction by abiotic or biotic stresses, VOC emissions were also studied. Both B. bassiana and T. afroharzianum induced a significant but transient (one to two-day-long) reduction of stomatal conductance, which may indicate rapid activation of defensive (rejection) responses, but also limited photosynthesis. At later stages, our results demonstrated a successful and complete plant colonization by B. bassiana, which induced higher photosynthesis and lower respiration rates, improved growth of roots, stems, leaves, earlier flowering, higher number of fruits and yield in tomato plants. Beauveria bassiana also helped tomato plants fight B. cinerea, whose symptoms in leaves were almost entirely relieved with respect to control plants. Less VOCs were emitted when plants were colonized by B. bassiana or infected by B. cinerea, alone or in combination, suggesting no activation of VOC-dependent defensive mechanisms in response to both fungi.

Immunity-associated volatile emissions of ß-ionone and nonanal propagate defence responses in neighbouring barley plants.

Plants activate biochemical responses to combat stress. (Hemi-)biotrophic pathogens are fended off by systemic acquired resistance (SAR), a primed state allowing plants to respond faster and more strongly upon subsequent infection. Here, we show that SAR-like defences in barley (Hordeum vulgare) are propagated between neighbouring plants, which respond with enhanced resistance to the volatile cues from infected senders. The emissions of the sender plants contained 15 volatile organic compounds (VOCs) associated with infection. Two of these, ß-ionone and nonanal, elicited resistance upon plant exposure. Whole-genome transcriptomics analysis confirmed that interplant propagation of defence in barley is established as a form of priming. Although gene expression changes were more pronounced after challenge infection of the receiver plants with Blumeria graminis f. sp. hordei, differential gene expression in response to the volatile cues of the sender plants included an induction of HISTONE DEACETYLASE 2 (HvHDA2) and priming of TETRATRICOPEPTIDE REPEAT-LIKE superfamily protein (HvTPL). Because HvHDA2 and HvTPL transcript accumulation was also enhanced by exposure of barley to ß-ionone and nonanal, our data identify both genes as possible defence/priming markers in barley. Our results suggest that VOCs and plant-plant interactions are relevant for possible crop protection strategies priming defence responses in barley.

European oak chemical diversity - from ecotypes to herbivore resistance.

Climate change is increasing insect pressure and forcing plants to adapt. Although chemotypic differentiation and phenotypic plasticity in spatially separated tree populations are known for decades, understanding their importance in herbivory resistance across forests remains challenging. We studied four oak forest stands in Germany using nontarget metabolomics, elemental analysis, and chemometrics and mapped the leaf metabolome of herbivore-resistant (T-) and herbivore-susceptible (S-) European oaks (Quercus robur) to Tortrix viridana, an oak pest that causes severe forest defoliation. Among the detected metabolites, we identified reliable metabolic biomarkers to distinguish S- and T-oak trees. Chemotypic differentiation resulted in metabolic shifts of primary and secondary leaf metabolism. Across forests, T-oaks allocate resources towards constitutive chemical defense enriched of polyphenolic compounds, e.g. the flavonoids kaempferol, kaempferol and quercetin glucosides, while S-oaks towards growth-promoting substances such as carbohydrates and amino-acid derivatives. This extensive work across natural forests shows that oaks' resistance and susceptibility to herbivory are linked to growth-defense trade-offs of leaf metabolism. The discovery of biomarkers and the developed predictive model pave the way to understand Quercus robur's susceptibility to herbivore attack and to support forest management, contributing to the preservation of oak forests in Europe.

Volatile organic compound patterns predict fungal trophic mode and lifestyle.

Fungi produce a wide variety of volatile organic compounds (VOCs), which play central roles in the initiation and regulation of fungal interactions. Here we introduce a global overview of fungal VOC patterns and chemical diversity across phylogenetic clades and trophic modes. The analysis is based on measurements of comprehensive VOC profiles of forty-three fungal species. Our data show that the VOC patterns can describe the phyla and the trophic mode of fungi. We show different levels of phenotypic integration (PI) for different chemical classes of VOCs within distinct functional guilds. Further computational analyses reveal that distinct VOC patterns can predict trophic modes, (non)symbiotic lifestyle, substrate-use and host-type of fungi. Thus, depending on trophic mode, either individual VOCs or more complex VOC patterns (i.e., chemical communication displays) may be ecologically important. Present results stress the ecological importance of VOCs and serve as prerequisite for more comprehensive VOCs-involving ecological studies.

Nitric oxide coordinates growth, development, and stress response via histone modification and gene expression.

Nitric oxide (NO) is a signaling molecule with multiple regulatory functions in plant physiology and stress response. In addition to direct effects on transcriptional machinery, NO executes its signaling function via epigenetic mechanisms. We report that light intensity-dependent changes in NO correspond to changes in global histone acetylation (H3, H3K9, and H3K9/K14) in Arabidopsis (Arabidopsis thaliana) wild-type leaves, and that this relationship depends on S-nitrosoglutathione reductase and histone deacetylase 6 (HDA6). The activity of HDA6 was sensitive to NO, demonstrating that NO participates in regulation of histone acetylation. Chromatin immunoprecipitation sequencing and RNA-seq analyses revealed that NO participates in the metabolic switch from growth and development to stress response. This coordinating function of NO might be particularly important in plant ability to adapt to a changing environment, and is therefore a promising foundation for mitigating the negative effects of climate change on plant productivity.

Heatwave frequency and seedling death alter stress-specific emissions of volatile organic compounds in Aleppo pine.

Biogenic volatile organic compounds (BVOC) play important roles in plant stress responses and can serve as stress indicators. While the impacts of gradual environmental changes on BVOCs have been studied extensively, insights in emission responses to repeated stress and recovery are widely absent. Therefore, we studied the dynamics of shoot gas exchange and BVOC emissions in Pinus halepensis seedlings during an induced moderate drought, two four-day-long heatwaves, and the combination of drought and heatwaves. We found clear stress-specific responses of BVOC emissions. Reductions in acetone emissions with declining soil water content and transpiration stood out as a clear drought indicator. All other measured BVOC emissions responded exponentially to rising temperatures during heat stress (maximum of 43 °C), but monoterpenes and methyl salicylate showed a reduced temperature sensitivity during the second heatwave. We found that these decreases in monoterpene emissions between heatwaves were not reflected by similar declines in their internal storage pools. Because stress intensity was extremely severe, most of the seedlings in the heat-drought treatment died at the end of the second heatwave (dark respiration ceased). Interestingly, BVOC emissions (methanol, monoterpenes, methyl salicylate, and acetaldehyde) differed between dying and surviving seedlings, already well before indications of a reduced vitality became visible in gas exchange dynamics. In summary, we could clearly show that the dynamics of BVOC emissions are sensitive to stress type, stress frequency, and stress severity. Moreover, we found indications that stress-induced seedling mortality was preceded by altered methanol, monoterpene, and acetaldehyde emission dynamics.

Protein expression plasticity contributes to heat and drought tolerance of date palm.

Climate change is increasing the frequency and intensity of warming and drought periods around the globe, currently representing a threat to many plant species. Understanding the resistance and resilience of plants to climate change is, therefore, urgently needed. As date palm (Phoenix dactylifera) evolved adaptation mechanisms to a xeric environment and can tolerate large diurnal and seasonal temperature fluctuations, we studied the protein expression changes in leaves, volatile organic compound emissions, and photosynthesis in response to variable growth temperatures and soil water deprivation. Plants were grown under controlled environmental conditions of simulated Saudi Arabian summer and winter climates challenged with drought stress. We show that date palm is able to counteract the harsh conditions of the Arabian Peninsula by adjusting the abundances of proteins related to the photosynthetic machinery, abiotic stress and secondary metabolism. Under summer climate and water deprivation, these adjustments included efficient protein expression response mediated by heat shock proteins and the antioxidant system to counteract reactive oxygen species formation. Proteins related to secondary metabolism were downregulated, except for the P. dactylifera isoprene synthase (PdIspS), which was strongly upregulated in response to summer climate and drought. This study reports, for the first time, the identification and functional characterization of the gene encoding for PdIspS, allowing future analysis of isoprene functions in date palm under extreme environments. Overall, the current study shows that reprogramming of the leaf protein profiles confers the date palm heat- and drought tolerance. We conclude that the protein plasticity of date palm is an important mechanism of molecular adaptation to environmental fluctuations.

Isoprene and ß-caryophyllene confer plant resistance via different plant internal signalling pathways.

Isoprene and other terpenoids are important biogenic volatile organic compounds in terms of atmospheric chemistry. Isoprene can aid plant performance under abiotic stresses, but the fundamental biological reasons for the high emissions are not completely understood. Here, we provide evidence of a previously unrecognized ecological function for isoprene and for the sesquiterpene, ß-caryophyllene. We show that isoprene and ß-caryophyllene act as core components of plant signalling networks, inducing resistance against microbial pathogens in neighbouring plants. We challenged Arabidopsis thaliana with Pseudomonas syringae, after exposure to pure volatile terpenoids or to volatile emissions of transformed poplar or Arabidopsis plants. The data suggest that isoprene induces a defence response in receiver plants that is similar to that elicited by monoterpenes and depended on salicylic acid (SA) signalling. In contrast, the sesquiterpene, ß-caryophyllene, induced resistance via jasmonic acid (JA)-signalling. The experiments in an open environment show that natural biological emissions are enough to induce resistance in neighbouring Arabidopsis. Our results show that both isoprene and ß-caryophyllene function as allelochemical components in complex plant signalling networks. Knowledge of this system may be used to boost plant immunity against microbial pathogens in various crop management schemes.

Ecotrons: Powerful and versatile ecosystem analysers for ecology, agronomy and environmental science.

Ecosystems integrity and services are threatened by anthropogenic global changes. Mitigating and adapting to these changes require knowledge of ecosystem functioning in the expected novel environments, informed in large part through experimentation and modelling. This paper describes 13 advanced controlled environment facilities for experimental ecosystem studies, herein termed ecotrons, open to the international community. Ecotrons enable simulation of a wide range of natural environmental conditions in replicated and independent experimental units while measuring various ecosystem processes. This capacity to realistically control ecosystem environments is used to emulate a variety of climatic scenarios and soil conditions, in natural sunlight or through broad-spectrum lighting. The use of large ecosystem samples, intact or reconstructed, minimizes border effects and increases biological and physical complexity. Measurements of concentrations of greenhouse trace gases as well as their net exchange between the ecosystem and the atmosphere are performed in most ecotrons, often quasi continuously. The flow of matter is often tracked with the use of stable isotope tracers of carbon and other elements. Equipment is available for measurements of soil water status as well as root and canopy growth. The experiments ran so far emphasize the diversity of the hosted research. Half of them concern global changes, often with a manipulation of more than one driver. About a quarter deal with the impact of biodiversity loss on ecosystem functioning and one quarter with ecosystem or plant physiology. We discuss how the methodology for environmental simulation and process measurements, especially in soil, can be improved and stress the need to establish stronger links with modelling in future projects. These developments will enable further improvements in mechanistic understanding and predictive capacity of ecotron research which will play, in complementarity with field experimentation and monitoring, a crucial role in exploring the ecosystem consequences of environmental changes.

Improving Air Quality by Nitric Oxide Consumption of Climate-Resilient Trees Suitable for Urban Greening.

Nitrogen oxides (NOx), mainly a mixture of nitric oxide (NO) and nitrogen dioxide (NO2), are formed by the reaction of nitrogen and oxygen compounds in the air as a result of combustion processes and traffic. Both deposit into leaves via stomata, which on the one hand benefits air quality and on the other hand provides an additional source of nitrogen for plants. In this study, we first determined the NO and NO2 specific deposition velocities based on projected leaf area (sV d) using a branch enclosure system. We studied four tree species that are regarded as suitable to be planted under predicted future urban climate conditions: Carpinus betulus, Fraxinus ornus, Fraxinus pennsylvanica and Ostrya carpinifolia. The NO and NO2 sVd were found similar in all tree species. Second, in order to confirm NO metabolization, we fumigated plants with 15NO and quantified the incorporation of 15N in leaf materials of these trees and four additional urban tree species (Celtis australis, Alnus spaethii, Alnus glutinosa, and Tilia henryana) under controlled environmental conditions. Based on these 15N-labeling experiments, A. glutinosa showed the most effective incorporation of 15NO. Third, we tried to elucidate the mechanism of metabolization. Therefore, we generated transgenic poplars overexpressing Arabidopsis thaliana phytoglobin 1 or 2. Phytoglobins are known to metabolize NO to nitrate in the presence of oxygen. The 15N uptake in phytoglobin-overexpressing poplars was significantly increased compared to wild-type trees, demonstrating that the NO uptake is enzymatically controlled besides stomatal dependence. In order to upscale the results and to investigate if a trade-off exists between air pollution removal and survival probability under future climate conditions, we have additionally carried out a modeling exercise of NO and NO2 deposition for the area of central Berlin. If the actually dominant deciduous tree species (Acer platanoides, Tilia cordata, Fagus sylvatica, Quercus robur) would be replaced by the species suggested for future conditions, the total annual NO and NO2 deposition in the modeled urban area would hardly change, indicating that the service of air pollution removal would not be degraded. These results may help selecting urban tree species in future greening programs.

Corrigendum: The Defense-Related Isoleucic Acid Differentially Accumulates in Arabidopsis Among Branched-Chain Amino Acid-Related 2-Hydroxy Carboxylic Acids.

[This corrects the article DOI: 10.3389/fpls.2018.00766.].

Metabolomic adjustments in the orchid mycorrhizal fungus Tulasnella calospora during symbiosis with Serapias vomeracea.

All orchids rely on mycorrhizal fungi for organic carbon, at least during early development. In fact, orchid seed germination leads to the formation of a protocorm, a heterotrophic postembryonic structure colonized by intracellular fungal coils, thought to be the site of nutrient transfer. The molecular mechanisms underlying mycorrhizal interactions and metabolic changes induced by this symbiosis in both partners remain mostly unknown. We studied plant-fungus interactions in the mycorrhizal association between the Mediterranean orchid Serapias vomeracea and the basidiomycete Tulasnella calospora using nontargeted metabolomics. Plant and fungal metabolomes obtained from symbiotic structures were compared with those obtained under asymbiotic conditions. Symbiosis induced substantial metabolomic alterations in both partners. In particular, structural and signaling lipid compounds increased markedly in the external fungal mycelium growing near the symbiotic protocorms, whereas chito-oligosaccharides were identified uniquely in symbiotic protocorms. This work represents the first description of metabolic changes occurring in orchid mycorrhiza. These results - combined with previous transcriptomic data - provide novel insights on the mechanisms underlying the orchid mycorrhizal association and open intriguing questions on the role of fungal lipids in this symbiosis.

Root isoprene formation alters lateral root development.

Isoprene is a C5 volatile organic compound, which can protect aboveground plant tissue from abiotic stress such as short-term high temperatures and accumulation of reactive oxygen species (ROS). Here, we uncover new roles for isoprene in the plant belowground tissues. By analysing Populus x canescens isoprene synthase (PcISPS) promoter reporter plants, we discovered PcISPS promoter activity in certain regions of the roots including the vascular tissue, the differentiation zone and the root cap. Treatment of roots with auxin or salt increased PcISPS promoter activity at these sites, especially in the developing lateral roots (LR). Transgenic, isoprene non-emitting poplar roots revealed an accumulation of O2- in the same root regions where PcISPS promoter activity was localized. Absence of isoprene emission, moreover, increased the formation of LRs. Inhibition of NAD(P)H oxidase activity suppressed LR development, suggesting the involvement of ROS in this process. The analysis of the fine root proteome revealed a constitutive shift in the amount of several redox balance, signalling and development related proteins, such as superoxide dismutase, various peroxidases and linoleate 9S-lipoxygenase, in isoprene non-emitting poplar roots. Together our results indicate for isoprene a ROS-related function, eventually co-regulating the plant-internal signalling network and development processes in root tissue.

The isoleucic acid triad: distinct impacts on plant defense, root growth, and formation of reactive oxygen species.

Isoleucic acid (ILA), a branched-chain amino acid-related 2-hydroxycarboxylic acid, occurs ubiquitously in plants. It enhances pathogen resistance and inhibits root growth of Arabidopsis. The salicylic acid (SA) glucosyltransferase UGT76B1 is able to conjugate ILA. Here, we investigate the role of ILA in planta in Arabidopsis and reveal a triad of distinct responses to this small molecule. ILA synergistically co-operates with SA to activate SA-responsive gene expression and resistance in a UGT76B1-dependent manner in agreement with the observed competitive ILA-dependent repression of SA glucosylation by UGT76B1. However, ILA also shows an SA-independent stress response. Nitroblue tetrazolium staining and pharmacological experiments indicate that ILA induces superoxide formation of the wild type and of an SA-deficient (NahG sid2) line. In contrast, the inhibitory effect of ILA on root growth is independent of both SA and superoxide induction. These effects of ILA are specific and distinct from its isomeric compound leucic acid and from the amino acid isoleucine. Leucic acid and isoleucine do not induce expression of defense marker genes or superoxide production, whereas both compounds inhibit root growth. All three responses to ILA are also observed in Brassica napus.

Sniffing fungi - phenotyping of volatile chemical diversity in Trichoderma species.

Volatile organic compounds (VOCs) play vital roles in the interaction of fungi with plants and other organisms. A systematic study of the global fungal VOC profiles is still lacking, though it is a prerequisite for elucidating the mechanisms of VOC-mediated interactions. Here we present a versatile system enabling a high-throughput screening of fungal VOCs under controlled temperature. In a proof-of-principle experiment, we characterized the volatile metabolic fingerprints of four Trichoderma spp. over a 48 h growth period. The developed platform allows automated and fast detection of VOCs from up to 14 simultaneously growing fungal cultures in real time. The comprehensive analysis of fungal odors is achieved by employing proton transfer reaction-time of flight-MS and GC-MS. The data-mining strategy based on multivariate data analysis and machine learning allows the volatile metabolic fingerprints to be uncovered. Our data revealed dynamic, development-dependent and extremely species-specific VOC profiles from the biocontrol genus Trichoderma. The two mass spectrometric approaches were highly complementary to each other, together revealing a novel, dynamic view to the fungal VOC release. This analytical system could be used for VOC-based chemotyping of diverse small organisms, or more generally, for any in vivo and in vitro real-time headspace analysis.

Origin of volatile organic compound emissions from subarctic tundra under global warming.

Warming occurs in the Arctic twice as fast as the global average, which in turn leads to a large enhancement in terpenoid emissions from vegetation. Volatile terpenoids are the main class of biogenic volatile organic compounds (VOCs) that play crucial roles in atmospheric chemistry and climate. However, the biochemical mechanisms behind the temperature-dependent increase in VOC emissions from subarctic ecosystems are largely unexplored. Using 13 CO2 -labeling, we studied the origin of VOCs and the carbon (C) allocation under global warming in the soil-plant-atmosphere system of contrasting subarctic heath tundra vegetation communities characterized by dwarf shrubs of the genera Salix or Betula. The projected temperature rise of the subarctic summer by 5°C was realistically simulated in sophisticated climate chambers. VOC emissions strongly depended on the plant species composition of the heath tundra. Warming caused increased VOC emissions and significant changes in the pattern of volatiles toward more reactive hydrocarbons. The 13 C was incorporated to varying degrees in different monoterpene and sesquiterpene isomers. We found that de novo monoterpene biosynthesis contributed to 40%-44% (Salix) and 60%-68% (Betula) of total monoterpene emissions under the current climate, and that warming increased the contribution to 50%-58% (Salix) and 87%-95% (Betula). Analyses of above- and belowground 12/13 C showed shifts of C allocation in the plant-soil systems and negative effects of warming on C sequestration by lowering net ecosystem exchange of CO2 and increasing C loss as VOCs. This comprehensive analysis provides the scientific basis for mechanistically understanding the processes controlling terpenoid emissions, required for modeling VOC emissions from terrestrial ecosystems and predicting the future chemistry of the arctic atmosphere. By changing the chemical composition and loads of VOCs into the atmosphere, the current data indicate that global warming in the Arctic may have implications for regional and global climate and for the delicate tundra ecosystems.